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The high parton density regime of the Quantum Chromodynamics (QCD), where 
the physics of parton saturation is expected to be dominant, is briefly discussed. Some 
phenomenological aspects of saturation are described, mainly focusing on possible sig- 
natures of the non-linear QCD dynamics in the heavy quark production in electron- 
proton/nucleus collisions. Implications of these effects in the heavy quark production in 
ultrapcriphcral heavy-ion collisions are also presented. 
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1. Introduction 

The understanding and analytic description of the QCD at high energies (small 
Bjorken x) has become an increasingly active subject of research, both from ex- 
perimental and theoretical points of view (For a recent review, see e.g. 0). These 
studies are mainly motivated by the violation of the unitarity (or Froissart-Martin 
bound El, which states that a to t < Cln 2 (s) at asymptotically large energies s) by 
the solutions of the linear perturbative DGLAP^and BFKL El evolution equations. 
Since these evolution equations predict that the cross section rises obeying a power 
law on energy, new dynamical effects associated with the unitarity restoration are 
expected to stop its further growth^. This expectation can be easily understood: 
while for large momentum transfer fej_, the BFKL equation predicts that the mech- 
anism g — > gg populates the transverse space with a large number of small size 
gluons per unit of rapidity (the transverse size of a gluon with momentum k^ is 
proportional to l/fej_), for small k^ the produced gluons overlap and fusion pro- 
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cesses, gg — > g, are equally important. Considering the latter process, the rise of 
the gluon distribution below a typical scale is slowed down, restoring the unitarity. 
That typical scale is energy dependent and is called saturation scale Q S at- The sat- 
uration momentum sets the critical transverse size for the unitarization of the cross 
sections. In other words, unitarity is restored by including non-linear corrections in 
the evolution equations*! Such effects are small for k± > Q 2 at and very strong for 



k 2 j_ < Qg at , leading to the saturation of the scattering amplitude. The magnitude 
of Qsat is associated to the behavior of the gluon distribution at high energies, and 
some estimates has been obtained. In general, the p redictions are Q sa t ~ 1 GeV 
at HERA/RHIC and Q sat - 2 - 3 GeV at LITC™ In particular, it has been ob- 
served that the ep deep inelastic scattering (DIS) data at low x can be successfully 
described with the help of the saturation model ^which incorporates the main 
characteristics of the high density QCD approaches I ' ' '. 

On the other hand, deep inelastic scattering on nuclei gives us a new possibility 
to reach high-density QCD phase without requiring extremely low values of x. The 
nucleus in this process serves as an amplifier for nonlinear phenomena. In order 
to understand this expectation and estimate the kinematic region where the high 
densities effects should be present, we can analyze the behavior of the function 
k(x,Q 2 ) = 3 2Q2 S X9A }tfP ^ i where xgA is the gluon distribution on the target A 
of transverse size Ra oc A 1 / 3 probed by a virtual probe of virtuality Q 2 . Such 
a function represents the probability of gluon-gluon interaction inside the parton 
cascade, and also is denoted the packing factor of partons in a parton cascade^. 
Considering that the condition k — 1 specifies the critical line, which separates 
between the linear (low parton density) regime k<1 and the high density regime 
k» 1, we can define the saturation momentum scale Q S at given by Q 2 &t (x\A) — 
37r 2 " 3 xgA ( x >Q^ x < A f> ; below which the gluon density reaches its maximum value 
(saturates). At any value of x there is a value of Q 2 — Q 2 &t (x) in which the gluonic 
density reaches a sufficiently high value that the number of partons stops to rise. 
This scale depends on the energy of the process [xg oc x~ x (A w 0.3)] and on the 
atomic number of the colliding nuclei [Ra a A 3 — > Q 2 oc As], with the saturation 
scale for nuclear targets larger than for nucleon ones. This result motivates more 
extensive studies of nuclear collisions and, in particular, of electron-nucleus collisions 
at high energies, where nuclear medium effects are reduced in comparison with AA 
collisions. 

In what follows we present a brief review of some signatures of the high parton 
density regime, with special emphasis on the heavy quark production in the sat- 
uration scenario. In the next section, we present a brief review of the saturation 
approaches for deep inelastic scattering process. In section [3] we discuss the prop- 
erty of geometric scaling in the inclusive charm production. Moreover, in Section 
01 we consider the heavy quark production in photonuclear process. In Section [5] 
the possibility of using ultraperipheral heavy ion collisions as a photonuclear col- 
lider is analyzed and some predictions for the saturation effects in the heavy quark 
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production are presented. Finally, in SectionElwe summarize our main conclusions. 



2. Overview on the Saturation Approaches 

We start from the space-time picture of the eletron-proton/nuclei processes The 
deep inelastic scattering ep(A) — > e + X is characterized by a large electron energy 
loss v (in the target rest frame) and an invariant momentum transfer q 2 = —Q 2 
between the incoming and outgoing electron such that x = Q 2 jlm^v is fixed (mjv 
is the target mass). In terms of Fock states we then view the ep(A) scattering as 
follows: the electron emits a photon (|e >^ \ej >) with E 1 = v and p\ 1 « Q 2 , 
after the photon splits into a qq Qe'y >— > \eqq >) and typically travels a distance 
l c w 1/m^x, referred as the coherence length, before interacting in the target. 
For small x, the photon converts to a quark pair at a large distance before its 
scattering. Consequently, the space-time picture of the DIS in the target rest frame 
can be viewed as the decay of the virtual photon at high energy into a quark- 
antiquark pair (color dipole), which subsequently interacts with the target. In the 
small x region, the color dipole crosses the target with fixed transverse distance r \ 
between the quarks. The interaction r )*p{A) is further factorized and is given byl^J 

al^ A) (x,Q 2 ) = J2 [ dzd 2 r ± \¥» T (z,r ± ,Q 2 )\ 2 4f(x,r x ), (1) 

where z is the longitudinal momentum fraction of the quark of flavor /. The photon 
wave functions ^l,t are determined from light cone perturbation theory (For a 
review see, e. g., Ref.Id). 

The dipole hadron (nucleus) cross section o~dip contains all information about 
the target and the strong interaction physics. Currently, the most complete QCD 
based effective theory that describes the physics of hadronic interactions at very 
high energies is the Color Glass Condensate (CGC) This represents the small x 
gluonic components in a hadronic wavefunction and is named so since: Color stands 
for the charge carried by the gluons; Glass stands for a clear separation of time 
scales between the fast and slow degrees of the wavefunction; Condensate stands for 
the high density of gluons which can reach values of order C(l/a s ). The regime of a 
CGC is characterized by the limitation on the maximum phase-space parton density 
that can be reached in the hadron/nuclear wavefunction (parton saturation) and 
very high values of the QCD field strength F^ v w l/a s . The large values of the gluon 
distribution at saturation suggest the use of semi-classical methods, which allow to 
describe the small- x gluons inside a fast moving nucleus by a classical color field 
EH. I n the CGC formalism a dip can be computed in the eikonal approximation, 
resulting 

o-dip{x, r x ) = 2 J d 2 b ± [1 - S (x, n, b±)] , (2) 

where S is the S-matrix element at fixed impact parameter b± which encodes all the 
information about the hadronic scattering, and thus about the non-linear and quan- 
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turn effects in the hadron wave function. It can be obtained by solving the functional 
evolution equation in the rapidity y = ln(l/a;) derive d by Jalilian-Marian, Iancu, 
McLerran, Weigert, Leonidov and Kovner (JIMWLK) . An equivalent equation 
was developed by Balitsky These equations form a set of coupled equations, 
and as such are very difficult to deal with analitically. An approximated equation 
which allows to deal with scattering at or near the unitarity limit was suggested by 
Kovchegov*2 and can be considered a mean field approximation for the Balitsky- 
JIMWLK equations. While the Kovchegov equation is not so complete it does have 
the advantage of being a precise nonlinear equation for a function. Many interest- 
ing limits of the Kovchegov equation has been understood by analytical methods, 
with the main properties for the S-matrix being: (a) for the interaction of a small 
dipole (r_i_ <C 1/Qsat), S(?-_l) « 1, which characterizes that this system is weakly 
interacting; (b) for a large dipole (r± ^> l/Q S at), the system is strongly absorbed 
which implies S(rjJ <C 1. This property is associate to the large density of saturated 
gluons in the hadron wave function. 

In our analysis we will consider the phenomenological saturation model proposed 
in Ref. ^ which encodes the main properties of the non-linear QCD approaches. In 
this model one has, 



CO 



= 1 — S (x, r±) ; S = exp 



(3) 



with a dip /co the scattering amplitude, averaged over all impact parameters 
and Q 2 at ~ A 2 e XXn i x °/ x ). The parameters of the model were constrained from the 
HERA small x data, coming out typical values of order 1-2 GeV 2 for the momentum 
scale. We have that when Q 2 at {x) r\ <C 1, the model reduces to color transparency, 
whereas as one approaches the region Q 2 at (x)r 2 L ss 1, the exponential takes care 
of resumming many gluon exchanges, in a Glauber-inspired way. Intuitively, this is 
what happens when the proton starts to look dark. Moreover, a smooth transition 
to the photoproduction limit is obtained with a modification of the Bjorken variable 
as x — ^ x —— (Q 2 + Am 2 )/W 2 p and the large x threshold corrections are accounted 
for by multiplying equation above by a factor (1 — x) n [n — 5(7) for light (heavy) 
quarks] . 

One of the shortcomings of the saturation model is that the impact parameters 
are averaged over. Thus there may be saturation at the center of the proton, but 
the averaging gives much weight to the edges of the proton, where saturation is 
not present. The implicit assumption in the approach is that the proton is treated 
as being homogeneous in the transverse plane. In such case, the impact parameter 
profile is given by the Heaviside function, s(bj_) = O (bo — fejj, and is considered to 
be peaked at central impact parameter, namely at b± = 0. Actually, this procedure 
is oversimplified and more realistic profiles can be considered. For phenomenological 
purposes a Gaussian or a hard sphere assumption are commonly taken into account. 
Recently, the impact parameter dipole saturation model 14 was developed, recover- 
ing the known Glauber-Mueller dipole cross section as well as the DGLAP evolution 
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has been included in the dipole cross section Using a simple parameterization 
for the S-matrix element in terms of the gluon distribution, those authors have 
obtained a smooth matching onto the DGLAP evolution, at least in the leading 
order DGLAP formalism, improving significantly the fit to the data in the large Q 2 
region. 

Therefore, despite the saturation model to be very successful in describing HERA 
data, its functional form is only an approximation of the theoretical non-linear QCD 
approaches. Recently, a lot of work has been done to found analytical solutions from 
the saturation formalisms. Along these lines, currently intense theoretical studies 
has been performed tow ards an understanding of the BFKL approach in the border 
of the saturation region ED l n partic ular , a parameterization for the dipole cross 
section has been implemented in Ref. ^( where this quantity was constructed to 
smoothly interpolate between the limiting behaviors analytically under control: the 
solution of the BFKL equation for small dipole sizes, r± <C l/<5sat(^), and the 
Levin- Tuchin law ^1 for larger ones, rj_ ^ l/Q S at(^)- The model has been used 
in phenomenological studies on vector meso n production 021 diffractive processes 
1^1 and longitudinal structure function .fJ^H at HERA, as well as neutrino-nucleon 
total cross section E2. 



3. Heavy quarks in lepton-hadron collisions - Geometric Scaling 

An important feature of the available saturation approaches is the prediction of 
the geometric scaling. Namely, the total "f*p cross section at large energies is not 
a function of the two independent variab les a : and Q 2 , but is rather a function of 
the single variable r = Q 2 /<9sat- I n Ref. the authors have demonstrated that 
the geometric scaling predicted at low momenta Q 2 < Q 2 at (x) is preserved by the 
BFKL evolution up to relatively large virtualities, within the kinematical window 
QLt < Q 2 < Qsat/ A QCD- As demonstrated in Ref. ESI the HERA data on the 
proton structure function F2 are consistent with scaling at x < 0.01 and Q 2 < 400 
GeV 2 . Similar behavior have been observed in exclusive processes I^H and in the 
nuclear case 12^1 

In this section we analyze the geometric scaling in the inclusive charm production 
From the experimental point of view, the HERA experiments have published 
data for the contribution of charmed meson production to the structure function 
i<2. This allows one to single out the charm contribution F£ to the total structure 
function and thus to investigate if the property of geometric scaling is also present in 
this observable. Before presenting our results, lets perform a qualitative analysis of 
the inclusive charm production using the saturation model El in order to shed light 
on the dipole configurations dominating the process in the relevant kinematical 
limits and show how the geometric scaling comes out. A characteristic feature in 
heavy quark production within the color dip ole approach is that the process is 
dominated by small size dipole configurations ^ . The overlap function weighting 
the dipole cross section is peaked at r± ~ l/m c ~ 0.1 fm even for sufficiently low 
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Fig. 1. Experimental data on inclusive charm production plotted versus the scaling variable. The 
curves are the saturation model (solid line) and symmetric saturation model (dashed lines). 



Q 2 values. As a consequence, charm production is dominated by color transparency 
and saturation effects are not important there, i.e. Odi P — Co Q 2 at (a;)T' 2 L /4. 

For the HERA kinematical region, we have Q 2 at ~ 1 GeV 2 , which implies that 
the relation Q 2 at < Q 2 + [i 2 is ever satisfied, where p? c = Am 2 . Consequently, we 
can define two kinematical regimes depending of the relation between Q 2 and /i 2 . 
For Q 2 3> [J. 2 we have scaling with logarithmic enhancement coming from aligned 
jet configurations, whereas for Q 2 <C /i 2 only symmetric dipole configurations con- 
tribute. Therefore, one obtains that the total cross section reads as, 

o& ~ ( 1 + ^ %) 0(Q 2 - Mc) + ^f^W 2 - Q 2 ) , (4) 

where the first term provides the behavior 1 /t at large r whereas the second term 
leads to a smooth transition down to the asymptotic (r-independent) behavior at 
small r. 

An analytical expression for the r dependence of the inclusive charm production 
can also be obtained in a less model de pen dent way. For this purpose we will make 
use of the symmetric saturation model 1^1, where the energy evolution of the pro- 
ton leads to the parton multiplication and the transverse momentum scale Q S a.t{x) 
appears. The main assumption is that the evolved proton can be described by a 
collection of independent dipoles at the time of the interaction whose sizes are dis- 
tributed around l/Q S at- The rate of growth of the parton densities is assumed to 
be Q 2 at (a;)/A 2 and the symmetry between low and high virtualities in jp interac- 
tions comes from the symmetry in the dipole-dipole cross section. In the HERA 
kinematical regime, i.e. Q 2 at < /i 2 , the inclusive charm production is given by, 



sat 



= A 2 ^ ^ _ 6Xp 



(l + log(r + T c ) 



(5) 



where r c = /i 2 /Q 2 at (a;) and the parameters N, A 2 and f(>,<) are taken from the 
data fit in Refill Here, we make the simplified assumption that the coupling with 
the dipole is flavor blind, in such way that N c£ — (2/5) A, with N being the global 
normalization describing Fi data. The factor 2/5 corresponds to the charge fraction 
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e c/ E e u d s + e c)- Once the parameters are fitted to proton structure function data, 
our prediction for the r dependence in the inclusive charm production is parameter 
free. 

In Fig. ^ we show the experimental data on the total cross section for the in- 
clusive charm production plotted versus scaling variable r, with Q sa t from the sat- 
uration model. We see the data exhibit geometric scaling for the whole Q 2 range, 
verifying a transition in the behavior on r of the cross section from a smooth de- 
pendence at small r and an approximated 1/r behavior at large r. The transition 
point is placed at fi 2 — 4m 2 , which takes values of order 10 GeV 2 for a charm 
mass m c = 1.5 GeV. This turns out in r ~ 10 since at HERA Q 2 at ~ 1 GeV 2 . 
The asymptotic 1/r dependence reflects the fact the charm production cross sec- 
tion scales as Q 2 at /Q 2 modulo a logarithmic correction ~ \n(Q 2 /fi 2 ), with energy 
dependence driven by the saturation scale. The mild dependence at r < /i 2 corre- 
sponds to the fact the cross section scales as Q 2 at /^c towards the photoproduction 
limit, but with the same energy behavior given by the saturation scale. In Fig. ^ 
we also found a symmetry between the regions of large and small r for the function 
y/r ofo t with respect the transformation r <-> 1/r in the whole region of r. The 
features present in the inclusive charm production data can be well reproduced in 
the phenomenological saturation model, corresponding to the solid curve in Fig. 
^ The symmetric saturation model also provides similar results, as shown in the 
dot-dashed lines. Disregarding the Glauber-like resummation in this model, the ex- 
pression gets simplified to a$£ t oc T _^ r [1 + log(r + r c )], and the symmetric pattern 
is easily verified. 

In the HERA kinematic domain the saturation momentum Q 2 at (a;) stays below 
the hard scale /i 2 = 4m 2 , implying that charm production probes mostly the color 
transparency regime and saturation corrections are not very important. However, 
as the saturation scale rises with the energy and the atomic number, we expect that 
at larger energies and nuclear collisions a new kinematic regime where Q 2 at (x) > /i 2 
will be probed. Recently, the production of open cha rm in heavy ion collisions in 
the CGC framework has been considered in Ref. 1221. The main prediction is the 
approximate scaling of the cross section with the number of participants (V part ) 
in the forward rapidity region, where the saturation scale exceeds the charm quark 
mass. This result is in contrast with usual expectation for a hard process of scaling 
with the number of collisions (iVcoii). These results provide a strong motivation for 
further investigations (See, e.g. Refs. I2U 31 | 32 |^ 

4. Heavy quark production in photonuclear process 

In this section, we report our investigations on the high energy hea vy quark pho- 
toproduction on nuclei targets using the saturation hypot hesi s 123 In particular, 
we study them considering the approach proposed in Ref. which extends the 
saturation model for scattering on nuclei and gives a reasonable parameter-free de- 
scription of the experimental data on nuclear structure function. In this model the 
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Fig. 2. The total nuclear cross section for charm and bottom photoproduction as a function of 
energy W 1 a f or distinct nuclei. 



dipole- nucleus cross section is given by^^, 

1 



oSf^fo r\,A)=2 / d 2 b ± U-exp 



--AT A (b x )a^ e °»(x, ri) 



(6) 



where 6j_ is the impact parameter of the center of the dipole relative to the center of 
the nucleus and the integrand gives the total dipole-nucleus cross section for a fixed 
impact parameter. The nuclear profile function is labeled by T^(b±). The above 
equation sums up all the multiple elastic rescattering diagrams of the qq pair and is 
justified for large coherence length, where the transverse separation r± of partons 
in the multiparton Fock state of the photon becomes as good a conserved quantity 
as the angular momentum, i. e. the size of the pair r± becomes eigenvalue of the 
scattering matrix. Here, the dipole cross section for the nucleon target (proton), 
u^fp n (x, is given by the saturation model [Eq. J2J]. The photoproduction 
cross section reads as, 

a to thA^QQX) = J dzd 2 r 1 _\^{z,r u Q 2 = 0)| 2 a^ lcus (x , r ± ) , (7) 

where the dipole-nucleus cross section is given by Eq. © and only the transverse 
wavefunction for the heavy quark- antiquark pair QQ contributes at Q 2 — > 0. 

In Fig. [21 are shown the results for the charm and bottom photoproduction cross 
section as a function of energy for different nuclei, including the proton case. The 
results present mild growth on W-yA at high energies stemming from the satura- 
tion model, whereas th e lo w energy region is consistently described through the 
large- a; threshold factor^!. For the proton, the experimental data from HERA and 
fixed target collisions are also included for comparison. The result for charm un- 
derestimates data by a factor 2 at W lp ~ 200 GeV, whereas is consistent with the 
measurements of bottom cross section. Concerning charm production, the measured 
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cross sections present a well known steeper behavior on energy even at electropro- 
duction, suggesting that further resummations in the original saturation model are 
needed in order to produce the larger growth on energy appearing in the charm 
measurements. We believe that the better result for bottom happens to be a mis- 
match between a large uncertainty in the experimental measurement and the lower 
bottom mass m& = 4.5 GeV considered here. For the nuclear case, we predict that 
their absolute values are rather large, reaching sa 2 • 10 3 and « 40 [ib for charm and 
bottom for lead at W^a = 10 3 GeV. 

The simple intuitive dipole approach considered above has a deep connection 
with a more general theoretical formalism. Namely, it is equivalent, at leading loga- 
rithmic approximation, to the high energy fc^-factorization (semihard) approach' 1 ^'. 
In this framework, the relevant QCD diagrams are considered with the virtualities 
and polarizations of the initial partons, carrying information on their transverse mo- 
menta. The scattering processes are described through the convolution of off-shell 
matrix elements with the unintegrated parton distribution, J r (x,kj_). The charac- 
teristic feature is the LO cross section in this approach resumming most of the NLO 
and even NNLO diagrams contributing to the process in the collinear formalism. As 
heavy quark photoproduction is concern ed, considering only the direct component 
of the photon, the cross section reads as l^^l 

vtotilA - QQX) = ^IS. J dz d i pi± S kl 



where D± = p 2 ± +mQ and D% = (kj_—p 1± ) 2 +m,Q. The transverse momenta of the 
heavy quark (antiquark) are denoted by p 1± and p 2± = (k± — Pi±), respectively. 
The heavy quark longitudinal momentum fraction is labeled by z. For the scale 
fj, in the strong coupling constant we use the prescription fi 2 = k\ + TOg. The 
unintegrated gluon distribution Fix, fej_) is directly related to the Fourier transform 
of the dipole-nucleon (nucleus) total cross section, as follows 

For the saturation model this quantity takes a simple analytical form, which 
reads for the proton and nucleus case as 1^1 

2n 2 a s \Qi A {x, b ± )J \ Q^ A [x, b±_) J 

where the expression for nuclei is obtained under the assumption of dominance of 
color transparency (small dipole configurations) in the dipole cross section, which 
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is correct for the particular case of heavy quark production. The nuclear satura- 
tion scale is given by Q 2 A (x, b±) = 4j4Ta(6x) fo Qsatfe)- We call attention to the 
scaling pattern on the variable r = k 2 j_/Q 2 A , which implies scaling on r also in 
the nuclear heavy quark production. In our further analysis, they will be compared 
with the simple ansatz J- nuc = dxGA ( x ,k±) ^ xG A (x,Q 2 ) being the nuclear 

gluon distribution (see f or details in the numerical calculations). 

On the other hand, in the collinear approach the cross section is given by a 
convolution between the partonic cross section for the subprocess jg — > QQ and 
the integrated gluon distr ibu tion for the nucleus xG A {x, Q 2 )- In our studies in Rcf. 

one considers the EKS^^and AG^I parameterizations for this distribution. The 
EKS parameterization was obtained from a global fit of the nuclear experimental 
data using the DGLAP evolution equations, which is a linear evolution equ atio n 
which does not consider dynamical saturation (high density) effects. In Ref. a 
procedure to include these effects in the nuclear gluon distribution was proposed, 
resulting in a paramerization for this distribution (AG parameterization), which 
also includes those present in the EKS parameterization. The AG parameterization 
predicts a stronger reduction of the growth of the gluon distribution at small values 
of x than the EKS one. 

When the nuclear photoproduction cross section of heavy quarks is computed 
considering these different approaches, we have obtained that the fcj_-factorization 
using EKS unintegrated gluon pdf (semihard approach) gives similar results to the 
collinear approach where nuclear effects (EKS parameterization) and high density 
corrections (AG parameterization) are taken into account. In particular, we have 
that the predictions using the AG parameterization in the collinear approach are 
similar to the semihard one, which docs not consider high density effects. This 
demonstrate that in this process we cannot distinguish if the modification in the 
behavior of the cross section is associated to high density effects in the collinear 
approach or a generalization of the factorization without high density effects in the 
unintegrated gluon distribution. On the other hand, if these effects are present and 
the factorization of the cross section is given by the fcj_-factorization, as is the case 
for the predictions from the saturation model, we have that the difference between 
the cross sections is large, which should allow to discriminate between the theo- 
retical approaches. Therefore, the nuclear cross section would provide a strong test 
concerning the robustness of the saturation approach in describing the observables. 
The situation is less clear comparing the semihard approach and the collinear one. 
One possible interpretation for this result is that the expected enhancement in the 
semihard approach, associated to the resummation of the (a s ln^^-) n in the co- 
efficient function, is not sizeable for inclusive quantities in the kinematic region of 
the future colliders. Probably, a more promising quantity to clarify this issue would 
be the transverse momentum distribution. In this case, the semihard approach 
seems to be in better agreement with experimental data in the pp collisions than 
the collinear approach 1*^1. 
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5. Heavy Quarks in Ultraperipheral Heavy Ion Collisions 

The studies of saturation effects in nuclear processes shown that future electron- 
nucleus colliders at HERA and RHIC, probably could determine whether parton 
distributions saturate and constrain the behavior of the nuclear gluon distribution 
in the full kinematical range. However, until these colliders become reality we need 
to consider alternative searches in the current and/or scheduled accelerators which 
allow us to constrain the QCD dynamics. Recently, we have analyzed the possibility 
of using ultraperipheral heavy ion collisions (UPC's) as a photonuclear collider. 
In particular, we have studied the heavy quark production I™ assuming distinct 
approaches for the QCD evolution. 

In heavy ion collisions the large number of photons coming from one of the 
colliding nuclei will allow to study photoproduction, with energies W^a reaching 
to almost 1 TeV for the LHC. The photonuclear cross sections are given by the 
convolution between the photon flux from one of the nuclei and the cross section for 
the scattering photon- nuclei, with the photon flux dN J^^ given by the Weizsacker- 
Williams method The final expression for the production of heavy quarks in 
ultraperipheral heavy ion collisions is given by, 

oo 

dN (cj) 

a AA^QQX 



NN ) = / du) —£ 1(7 ^QQx{^ A = 2u>^S^) (10) 



where to is the c.m.s. photon energy, to m i n — M'\ Tt jA r )i 1 m, v and v / Snn is the c.m.s 

WW 

energy of the nucleus- nucleus system. The Lorentz factor for LHC is 7l = 2930, 
giving the maximum c.m.s. jN energy W~/A ~ 950 GeV. The requirement that 
photoproduction is not accompanied by hadronic interaction (ultraperipheral colli- 
sion) can be done by restricting the impact parameter b to be larger than twice the 
nuclear radius, Ra — 1.2 A 1 / 3 fm. An analytic approximation for AA collisions can 
be obtained using as integration limit b > 2Ra, producing 



dN(uj) 2 Z 2 a em \_ ,_s r _s rf ,„ 2 ,_s ^ 2 /-n\ 
- ——— V K (v) Ki (rj) + — (K 1 (77) - K (77)) 



(11) 



where fj = 2ojRa/jl and Kq^(x) are the modified Bessel functions. The typical 
values of x which will be probed in ultraperipheral heavy ion collisions, are given by 
x = (AlQQ/2p)e~ y , where Mqq is the invariant mass of the photon-gluon system 
and y the center of momentum rapidity. For Pb + Pb collisions at LHC energies 
the nucleon momentum is equal to p = 2750 GeV; hence x = (A/qq/5500 GeV)e~ y . 
Therefore, the region of small mass and large rapidities probes directly the high 
energy (small x) behavior of the QCD dynamics present in the 7 A cross section. 
This demonstrates that ultraperipheral heavy ion collisions at LHC represents a 
very good tool to constrain the high energy regime of the QCD dynamics. 

In what follo ws w e summarize our analysis on photonuclear production of heavy 
quarks at UPC's In order to do this, one considers the available high energy 
approaches. Namely, the usual collinear approach, the semihard formalism and the 
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Fig. 3. Rapidity distribution for charm production in ultraperipheral heavy ion collisions. 

phenomcnological saturation reviewed in the previous section. As an additional 
analysis, we also consider the Color Glass Condensate (CGC) formalism. This puts 
into perspective the most recent representative high energy approaches and allow 
us to find out the observables and/or distributions which could disentangle them at 
the planne d co lliders. 

In Refs.EDIthe heavy quark production in UPC's has been analyzed in the CGC 
formalism. In Ref. GUI we have improved that analysis using a realistic photon flux 
and a color field correlator including quantum radiation effects. The differential 
cross section on rapidity reads as GUI 

±A^99E ^ u djpa^ + r # d{k±) L + 4 ( fe x ~f>) arcth *x 1 f 
dY dui 2 it J I k\ u% u% I 



where we define the rapidity Y = hx(l/x) = ln(2 lu 7 l / Auiq) ■ The quark charge is 

labeled as cq and one uses the notation p,Q = \Jk]_+ 4niQ. In Ref. ^| we obtained 
the following analytical expression for the color field correlator, considering that it 
is directly related to the Fourier transform of the dipole-nucleus total cross section, 

4tt \ / k 2 . 



c ^^{m^)^\-m^)- (12) 

where we have assumed Q 2 A (x) = A 1 / 3 Q 2 at (x). We believe that this input is more 
suitable for realistic computations because it includes quantum evolution in the 
formalism and reproduces most part of the phenomenological features of the sat- 
uration model for the nucleus case. It is worth mentioning the direct relation be- 
tween the color field correlator and the unintegrated gluon distribution, given by 
T(x, k x ) = (3i?y87r 2 a s ) k\ C (x, fe±). 

In Fig. 13 is shown the charm rapidity distribution for the distinct high en- 
ergy approaches considered. The collinear result is denoted by the long-dashed 
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curves, where we have employed the EKS98 parameterization for the collinear nu- 
clear gluon function (See also Ref . I^lf) . The solid and dotted lines label the semi- 
hard (/c^-factorization) results, where one has used the unintegrated gluon distri- 
bution J- nuc (x, k\; A) discussed in the previous section. Two possibilities for the 
nucleon gluon distribution were considered: (I) GRV94(LO) - solid line - and (II) 
GRV98(LO) - dotted line. The saturation model results are denoted by the dot- 
dashed line. The CGC prediction is denoted by the dashed line. We have that the 
predictions for the collinear approach and the semihard formalism are similar and 
give somewhat larger values than the saturation and CGC results. One possible 
interpretation for the similarity between the predictions of the semihard approach 
and the collinear one is that the expected enhancement in the fcj_-factorization for- 
malism i s no t sizeable for inclusive quantities in the kinematic region of the future 
colliders ^1 Our phenomenological ansatz within the CGC formalism gives similar 
results as the saturation model, but should be noticed that the physical assump- 
tions in those models are distinct. While the saturation model considers multiple 
scattering on single nuclcons, our expression for the dipole-nucleus cross section in 
the CGC formalism assumes scattering on a black area filled by partons coming 
from many nucleons. 

Let us present the numerical calculation of their total cross section at UPC's. 
We focus mostly on LHC domain where small values of x would be probed. The 
results are presented in Tabled The collinear approach gives a larger rate, followed 
by the semihard approach. The saturation model and CGC formalisms give similar 
results, including a closer ratio for charm to bottom production. Concerning the 
CGC approach, our phenomenological educated guess for the color field correlator 
seems to produce quite reliable estimates. Therefore, the photonuclcar production 
of heavy quarks allow us to constraint already in the current nuclear accelerators 
the QCD dynamics since the main features from photon-nuclei collisions hold in the 
UPC reactions. Our results shown that an experimental analysis of this process can 
be useful to constrain the QCD dynamics at high energies. 

6. Summary 

The perturbative QCD has furnished a remarkably successful framework to inter- 
pret a wide range of high energy lepton-lepton, lepton-hadron and hadron-hadron 
processes. Through global analysis of these processes, detailed information on the 
parton structure of hadrons, especially the nucleon, has been obtained. The existing 
global analysis have been performed using the standard DGLAP evolution equa- 
tions. However, in the small x region the DGLAP evolution equations are expected 
to breakdown, since new dynamical effects associated to the high parton density 
must occur in this kinematical region. 

Research in the field of QCD at high parton density deals both with fundamental 
theoretical issues, such as unitarity of strong interactions at high energies, and with 
the challenge of describing experimental data coming, at present, from HERA and 
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Table 1. The photonuclear heavy quark total cross sections for UPC's at 
LHC. 



QQ 


Collinear 


SAT-MOD 


SEMIHARD I (II) 


CGC 


cc 


2056 mb 


862 mb 


2079 (1679.3) mb 


633 mb 


bb 


20.1 mb 


10.75 mb 


18 (15.5) mb 


8.9 mb 



RHIC and expected exciting physics of forthcoming experiments at LHC. Over the 
past few years much theoretical effort has been devoted towards the understanding 
of the growth of the total scattering cross sections with energy. These studies are 
mainly motivated by the violation of the unitarity (or Froissart) bound by the 
solutions of the linear perturbative DGLAP and BFKL evolution equations. Since 
these evolution equations predict that the cross section rises obeying a power law 
of the energy, violating the Froissart bound, unitarity corrections are expected to 
stop its further growth. 

In this paper we have presented a brief review of the basic concepts present in 
the high density approaches and discussed some aspects of the rapidly developing 
field of QCD at high parton density in ep, eA and AA collisions. The successful 
description of all inclusive and diffractive deep inelastic data at the collider HERA, 
as well as some recent results from RHIC, by saturation models suggests that these 
effects might become important in the energy regime probed by current colliders. 
In particular, the remarkable property of geometric scaling verified in the data 
indicate that the experiments are in a kinematical region which probe QCD in the 
non-linear regime of high parton density. These results show that the transition 
between the linear and non-linear regimes in eA processes at high energies will 
occur in a perturbative regime, justifying perturbative QCD approaches. Our recent 
studies shown that an alternative for eA colliders is the study of saturation effects 
in ultraperipheral heavy ion collisions. 
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